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Clusters of galaxies generally form by the gravitational merger of smaller clusters and groups. 
Major cluster mergers are the most energetic events in the Universe since the Big Bang. 
Mergers drive shocks into the intracluster gas, and these shocks heat the intracluster gas, 
and should also accelerate nonthermal relativistic particles. The X-ray signatures of the 
thermal effects of merger shocks will be discussed. X-ray observations of shocks can be used 
to determine the geometry and kinematics of the merger. As a result of particle acceleration in 
shocks, clusters of galaxies should contain very large populations of relativistic electrons and 
ions. Electrons with Lorentz factors 7 ^ 300 (energies E = 'ym^c^ ~ 150 MeV) are expected 
to be particularly common. Observations and models for the radio, extreme ultraviolet, hard 
X-ray, and gamma-ray emission from nonthermal particles accelerated in these shocks will also 
be described. The predicted gamma-ray fluxes of clusters should make them easily observable 
with GLAST. 



1 Introduction 



Major cluster mergers are the most energetic events in the Universe since the Big Bang. Cluster 
mergers are the mechanism by which clusters are assembled. In these mergers, the subclusters 
collide at velocities of ~2000 km/s, and shocks are driven into the intracluster medium. In major 
mergers, these hydrodynamical shocks dissipate energies of ~ 3 x 10^^ ergs; such shocks are the 
major heating source for the X-ray emitting intracluster medium. The shock velocities in merger 
shocks are similar to those in supernova remnants in our Galaxy, and we expect them to produce 
similar effects. Mergers shocks should heat and compress the X-ray emitting intracluster gas, 
and increase its entropy. We also expect that particle acceleration by these shocks will produce 
nonthermal electrons and ions, and these can produce synchrotron radio, inverse Compton (IC) 
EUV and hard X-ray, and gamma-ray emission. 

Hydrodynarnical simulations of cluster formation and evolution have shown the importance 
of merger shocksE^'Ej The evolution of the structure of merger shocks is illustrated in Figure fflP 




Figure 1: The results of a hydro simulation of a symmetric, off-center merger by Richer and SarazinEH The colors 
show the temperature, while the contours are the X-ray surface brightness. Initially, the shocked region is located 
between the two subcluster centers. Later, the main merger shocks propogate to the outer parts of the cluster, and 
other weaker shocks also occur. By the end of the simulation, the cluster is beginning to return to equilibrium. 

which shows an off-center merger between two symmetric subclusters. At early stages in the 
merger (prior to the first panel), the shocked region is located between the two subcluster centers 
and is bounded on either side by two shocks. At this time, the subcluster centers, which may 
contain cooling cores and central radio sources, are not affected. Later, these shocks sweep over 
the subcluster centers (the first panel); the siiiYJijal of central cool cores depends on how sharply 
peaked the potentials of the subclusters arepiJE3 The main merger shocks pass into the outer 
parts of the merging system (panel 2), and secondary shocks may appear in the inner regions 
(panel 3). Eventually, the cluster begins to return to equilibrium (panel 4). 



2 Thermal Effects of Merger Shocks 



2.1 Shock Kinematics 

Merger shocks heat and compress the intracluster gas, and these effects can be used to determine 
the geometry and kinematics of the merger. ASCA X-ray temperature maps and ROSAT images 
have been used in an initial effort to apply this techniqueEll The cluster containing the radio 
source Cygnus A is a particularly simple case (Fig. |2|). This appears to be a fairly symmetric 
merger with a low impact parameter. The merger is at an early phase, with the merger shocks 
being located between the two subcluster centers. A hydro/N-body simulation of the mergeiEll 
is shown at the right in Fig. ^ (not to the same scale). Presumably, the fact that the merger 
shocks have not yet passed through the subcluster centers is the reason why the merger hasn't 
disrupted the Cygnus A radio source or the surrounding cooling flow. 

The simple geometry of this merger makes it easy to apply the shock jump conditions to 
determine the merger velocity. From the Rankine-Hugoniot jump conditions, the velocity change 




Figure 2; (a) ROSAT PSPC contours of the Cygnus A cluster. The radio galaxy is located at the center of 
the cooling flow peak in the southeast subcluster. A second subcluster is centered on the X-ray peak 11' to the 
northwest (labeled 6). The color shows the ASCA temperatures (red bpkig hot, blue being cool), and indicates 
that a merger shock regioftiis located between the two subcluster centersUil (b) A hydro/N-body simulation of the 
Cygnus A cluster mergerpZI not to the same scale as panel (a) but rotated to the same orientation. 



across the merger shock is 
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where Tq and Ti are the pre- and post-shock temperature, /i is the mean mass per particle, and 
the shock compression r is given by 
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For a symmetric merger, the merger velocity of the two subclusters is just Aud = 2Aush- When 
the ASCA temperatures are used, this gives Aud ~ 2200 km/s. The radial velocity distribution 
of the galaxies in this cluster is bimodall^ and consistent with a merger velocity of 2400 km/s. 
Interestingly, the collision velocity found above is close to the free-fall velocity of --^2200 
km/s that the two subclusters should have achieved had they fallen from a large distance to their 
observed separation. This consistency suggests that the shock energy is effectively thermalized, 
and that a major fraction does not go into turbulence, magnetic fields, or cosmic rays. 



2.2 Nonequilihrium Effects 

Cluster merger are expected to produce collisionless shocks, as occurs in supernova remnants. 
As such, nonequilibrium effects are expected, including nonequipartition of electrons and ions 
and nonequilibrium ionizationEll'E3 In fact, observations of a central shock in the spectacular 
merging cluster Abell 3667 show an apparent lag of ~200 kpc between the shock position and 
the peak in the electron temperature tiJ This is consistent with weak electron heating at the 
shock, and the timescale for Coulomb heating. However, recent Chandra observations of this 
cluster suggest a more complex situationcj 
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Figum 3: The Chandra X-ray image of the central region of the merging cluster Abell 2142 from Markevitch 
et alfc3 Note the two bow-shaped discontinuities to the northwest and just south of the brightest region of the 
cluster. 

2.3 Early Chandra Results 

Although ASCA and ROSAT observations provide some evidence for merger shocks, they lack 
the combination of spatial and spectral resolution over a hard X-ray band needed to really 
determine shock structures accurately. Chandra has superb spatial resolution, while XMM has 
a wider field of view and a very large collecting area. Astro-E would have provided very high 
spectral resolution, sufficient to directly detect the Doppler motions in the gas. Hopefully, an 
equivalent capability will be provided in the not too distant future. 

A number of merging clusters were included in the early observations with Chandra taken for 
calibration and science verification purposes, and several of these have shown sharp structures 
which initially appeared to be shocks. An example is Abell 2142 (Fig. ^, for which a careful 
analysis has been given by Markevitch et alo Even the raw Chandra image shows two bow- 
shaped structures, to the south and northwest of the center. However, the Chandra spectra show 
that these are not shocks. As one moves inward and the gas density increases abruptly, the gas 
temperature decreases in such a way that the pressure is fairly continuous. Thus, the specific 
entropy is actually lower in the denser "shocked" region than in the lower density "pre-shock" 
region. These are not shocks, but rather contact discontinuities between higher density cool gas 
and lower density hot gas. Markevitch et alE3 argue that the merger shocks have already passed 
though the lower density gas in these regions, and that the denser cooler gas regions are the 
remnants of the cooling cores of the original subclusters, which are now plowinffthrough the 
lower density shocked gas. A very similar situation has been found in Abell 3667^3 



3 Particle Acceleration and Nonthermal Emission 

3.1 Shock Acceleration 

Radio observations of supernova remnants indicate that shocks with v ^ 10^ k m/s convert at 
least a few percent of the shock energy into the acceleration of relativistic electronsB^ Even more 
energy may go into relativistic ions. While the merger shocks in cluster have lower Mach numbers 
and compressions than the blast waves in young supernova remnants, it is probable that shock 



Figure 4: The loss time scale for relativistic electrons in a cluster of galaxies under typical intracluster conditions 
as a function of their Lorentz factor 7. The particles have energies of -E = ^m^c^ . 

acceleration operates efficiently in clusters as well. Given that all of the thermal energy of the 
intracluster gas in clusters is due to shocks with such velocities, it seems likely that relativistic 
electrons with a total energy of ^ lO^'^ ergs are produced in clusters, with perhaps even higher 
energies in ions. 

Clusters are also very good storage locations for cosmic rays. Under reasonable assumptions 
for the diffusion coefficient, particles with energies ^10^ GeV have diffusion times which are 
longer than the Hubble timeBEJ Figure ^ gives the loss time scales for relativistic electrons 
under typical intracluster conditions. Although high energy electrons lose energy rapidly due to 
IC and synchrotron emission, electrons with Lorentz factors of 7 ~ 300 (energies ~ 150 MeV) 
have long lifetimes which approach the Hubble timeulrcj Thus, clusters of galaxies can retain low 
energy electrons (7 ~ 300) and nearly all cosmic ray ions for a significant fraction of a Hubble 
time. 

Recently, I have calculated models for the rel ativ istic electrons in clusters, assuming they are 
primary electrons accelerated in merger shocksrHrn An alternative theory is that these particles 
are secondaries produced by the decay of charged mesons generated in cosmic ray ion collisions Li 
Two recent cluster merger simulations have included particle acceleration approximatelyS'c; 
Their conclusions are very similar to mine based on simpler models. The populations of cosmic 
ray electrons in clusters depends on their merger histories. Because low energy electrons have 
long lifetimes, one expects to find a large population of them in most clusters (any cluster which 
has had a significant merger since z ~ 1). On the other hand, higher energy electrons [E '^1 
GeV) have short lifetimes (shorter than the time for a merger shock to cross a cluster). Thus, 
one only expects to find large numbers of higher energy electrons in clusters which are having 
or have just had a merger. 

Fig. |5|(a) shows the electron spectrum in a cluster with a typical history. Most of the 
electron energy is in electrons with 7 ~ 300, which have the longest lifetimes. These electrons 
are produced by mergers over the entire history of the cluster. This cluster also has a small 
ongoing merger which produces the high energy tail on the electron distribution. 

Most of the emission from these electrons is due to IC, and the resulting spectrum is shown in 
Fig. |^(b). For comparison, thermal bremsstrahlung with a typical rich cluster temperature and 
luminosity is shown as a dashed curve. Fig. |5|(b) shows that clusters should be strong sources of 
extreme ultraviolet (EUV) radiation. Since this emission is due to electrons with^ ^ 300 which 
have very long lifetimes, EUV radiation should be a common feature of clusterscfl 

In clusters with an ongoing merger, the higher energy electrons will produce a hard X-ray 
tail via IC scattering of the Cosmic Microwave Background (CMB); the same electrons will 
produce diffuse radio synchrotron emission. 
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Figure 5: (a) A typical model for the relativistic electron population in a cluster of galaxies. The lower energy 
electrons are due to all of the mergers in the cluster history, while the high energy electrons are due to a small 
current merger, (b) The IC spectrum from the same model (solid curve). The dashed curve is a 7 keV thermal 
bremsstrahlung spectrum. 



3.2 EUV/Soft X-ray Emission 

Excess EUV emission has apparently been detected with the EUVE sat ellite in six clusters 
(Virgo, Coma, Abeh 1795, Abell 2199, Abell 4038, & Abell 4059) BSSiliiiiiiaiii In fact, 
the EUVE satellite appears to have detected all of the clusters it observed which are nearby, 
which have long integration times, and which lie in directions of low Galactic column where 
detection is possible at these energies. However, the EUV detections and claimed properties of 
the clusters remain quite controversial p'tfa'H 

The EUV observations suggest that rich clusters generally have EUV luminosities of ~10'*^ 
ergs/s, and have spectra which decline rapidly in going from the EUV to the X-ray band. The 
EUV emission is very spatially extended. The last point is illustrated in Figure 0. The data 
points show the ratio of the EUV to thermal X-ray emission in the cluster Abell 1795. The ratio 
increases rapidly with increasing radius. 

While it is possible that the EUV emission may be thermal in originpja I believe that it is 
more likely that this emission is due to inverse Compton scattering (IC) of CMB photons by 
low energy relativistic electronsO'&Ejo In this model, the EUV would be produced by electrons 
with energies of ~150 MeV (7 ~ 300; Fig. |5|). As noted above, these electrons have lifetimes 
which are comparable to the Hubble time, and should be present in essentially all clusters. This 
can explain why EUV emission is common. To produce the EUV luminosities observed, one 
needs a population of such electrons with a total energy of ~10^^ ergs, which is about 3% of the 
typical thermal energy content of clusters. This is a reasonable acceleration efficiency for these 
particles, given that both the thermal energy in the intracluster gas and the relativistic particles 
result from merger shocks. The steep spectrum in going from EUV to X-ray bands is predicted 



by this model [Fig. |5Kb)], it results from the rapid increase in losses (oc 7 ) for particles as the 
energy increases above 7 ~ 300 (Fig. ^). 

The broad spatial distribution of the EUV is also naturally explained by the density depen- 
dence of IC emission. The thermal emission which produces the bulk of the X-ray emission in 
clusters is due to collisions between thermal electrons and ions; thus, it declines with the square 
of the density as the radius increases. On the other hand, IC emission is due to collisions between 
cosmic ray electrons and CMB photons; since the CMB energy density is extremely uniform, IC 
EUV emission varies with a single power of density, rather than density squared. This simple 
difference can explain why the EUV is more extended than the thermal X-ray emission. As 
a example, the solid curve in Figure ^ shows the predicted ratio of EUV to X-ray emission in 
Abell 1795, if the pressure in cosmic ray electrons varies in proportion with the thermal gas 




radius (arcmm) 

Figure 6: The data points give the observecj-^atio of the EUV to thermal X-ray surface brightness versus radius 
in the cluster Abell 1795 from Mittaz et alLJ The EUV data is from,EUVE, while the X-ray data is from the 
ROSAT PSPC. The solid curve is the IC model from Sarazin & Lieufj assuming that the pressure in cosmic ray 
electrons in this cluster varies in proportion with the thermal pressure in the X-ray emitting gas, as determined 
with the ROSAT PSPC. 

pressure, and the latter is determined from the X-ray observations E3 

3.3 Hard X-ray Tails 

If clusters contain higher energy relativistic electrons with 7 ~ 10^, these particles will produce 
hard X-ray emission by IC scattering, and radio synchrotron emission depending on the intra- 
cluster magnetic field. Since these higher energy electrons have short lifetimes, they should only 
be present in clusters with evidence for a recent or ongoing merger. 

Hard X-ray emission in excess of the thermal emission and detected as a nonthermal tail at 
energies ;^20 keV has been seen in at least two clusters. The Coma cluster, which is undergoing at 
least one merger and which has a radio halo, was detected with both BeppoSAX and RXTEGJ'E3 
BeppoSAX also has detected Abell 2256^^ another merger cluster with strong diffuse radio 
emission. BeppoSAX may have detected Abell 2199p3 although I believe the evidence is less 
compelling for this case. A nonthermal hard X-ray detection of Abell 2199 would be surprising, as 
this cluster is very relaxed and has no radio halo or relicBa An alternative explanation of the hard 
X-ray tails is that they might be due to nonthermal bremsstrahlungpj which is bremsstrahlung 
from nonthermal electrons with energies of 10-1000 keV which are being accelerated to higher 
energies. 

3.4 Radio Halos and Relics 

A number of clusters of galaxies are known to contain large-scale diffuse radio sources which have 
no obvious connection to individual galaxies in the cluster£3 These sources are referred to as radio 
halos when they appear projected on the center of the cluster, and are called relics when they are 
found on the cluster periphery. The best known and studied radio halo is in the Coma clusters 
These radio halos and relics are relatively uncommon; about 40 are known at the present time. 
They are diffuse radio sources with very steeply declining radio spectra. In all cases of which 
I am aware they have been found in clusters which show significant evidence for an ungoing 
mergerEHIOllj This suggests that the relativistic electrons are accelerated by merger shocks. 
In the early stages of mergers, halos are often found on the border between the subclusters, 
where the cluster gas is first being shocked (e.g., Abell 85) e3 In more advanced mergers, more 
conventional centrally located halos (e.g.. Coma) and peripheral relics (e.g., Abell 3667) are 
found. Abell 3667 provides a spectacular case of two bow-shaped radio relics at opposite sides 
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Figure 7: The predicted gamma-ray spectrum for the Coma cluster, including electron bremsstrahlung and 7r° 
decay from ionsu3 

of a merging clusterpl and located at the positions where merger shocks are predictedESI 

Although radio halos and relics have been studied for some time, lately the number of 
known cases has increased dramatically. This is largely due to the radio surveys (like NVSS and 
WENSS) and new low frequency radio instruments. Perhaps the largest number of new sources 
have come from the NVSSsurvey and VLA pointed observation follow up effort by Giovannini, 
Feretti, and collaboratorstHI Other new detections have involved other VLA surveyspl the VLA 
at 74 MHz@ SUMSS and MOST@ the ATCA@ and WENSsS 

3.5 Predicted Gamma-Ray Em,ission 

Relativislic-^jelectrons and ions in clusters are also expected to produce strong gamma-ray 
emission tll'BQ'E^ The region near 100 MeV is particularly interesting, as this region includes 
bremsstrahlung from the most common electrons with 7 ~ 300, and gamma-rays from ions 
produced by tt° decay. Both of these processes involve collisions between relativistic particles 
(electrons for bremsstrahlung, ions for tt° emission) and thermal particles, so they should both 
vary in the same way with density in the cluster. Thus, the ratio of these two spectrally distin- 
guishable emission processes should tell us the ratio of cosmic ray ions to electrons in clustersPE3 
Figure ^ shows the predicted gamma-ray spectrum for the Conm cluster, based on a model 
which reproduces the observed EUV, hard X-ray, and radio emissipnc3 The observed upper limit 
from CGO/EGRET is <4 x 10"^ cts/cmVs for E > 100 MeVP while the predicted value for 
this model is ~2 x 10^^ cts/cm^/s. The EGRET upper limit already shows that the ratio of 
ions to electrons cannot be too large (;^30)pe3) The predicted fluxes are such that many nearby 
clusters should be easily detectable with GLAST. 
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